A two-level atom coupled to a controllable squeezed vacuum field reservoir 
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The dissipative and decoherence properties of the two-level atom interacting with the squeezed 
vacuum field reservoir are investigated based on the nonautonomous master equation of the atomic 
density matrix in the framework of algebraic dynamics. The nonautonomous master equation is 
converted into a Schrodinger-like equations and its dynamical symmetry is found based on the left 
and right representations of the relevant algebra. The time-dependent solution and the steady 
solution are obtained analytically. The asymptotic behavior of the solution is examined and the 
approach to the equilibrium state is proved. Based on the analytic solution the response of the 
, system to the squeezed vacuum field reservoir is studied numerically. 

^ 1 PACS numbers: 03.65.Fd, 03.65.Yz, 42.50.Lc 

I. INTRODUCTION 

The fundamental property of the squeezed state is that the quantum fluctuations in one quadrature component 
of the field can be reduced heavily. Followed the early works Q , much attention has been given to this specific 
state. The first experimental result for the generation of the squeezed state was reported by Slusher et al. 3] with the 
scheme of 4-wave mixing in an optical cavity. Currently a successful scheme for generating squeezed light can also be 
based on a parametric oscillator or parametric down converter ^[. Recently, due to its potential applications in the 
fields of quantum measurement, optical communication, and quantum information, the squeezed vacuum state has 
been extensively studied [H, ||| . A natural problem is what effect on physical systems can be induced by the squeezed 
vacuum. The squeezed light field will generally be characterized as a non-stationary reservoir which contains phase 
• dependent features in the correlation function between pairs of photons. When the bandwidths of the squeezed lights 
are not too small, they can be treated as Markovian reservoirs, and the master equations of the reduced system can 
(-h , be obtained based on Markovian approximation. 

Master equations are of fundamental importance in the treatment of dissipation and decoherence of open quantum 
systems. The common feature of the quantum master equations is the existence of the sandwich terms of the Liouville 
operators where the reduced density matrix of the system is in between some quantum excitation and de-excitation 
operators. So it is very difficult to get the exact analytical solution of the master equation, only simple cases such as 
a single mode of cavity field coupled to a vacuum state reservoir (T — 0) or stationary regime properties are the ones 
analytically treated Q. fnstead, the master equations are normally converted into c-number equations in the coherent 
state representation-the Fokker-Planck equation 0,0- On the other hand, with the development of the so called 
quantum engineering, the man-made nonautonomous quantum systems where the system parameters set by people 
for controlling of the system are time-dependent, become more and more important, ft would be very desirable to get 
the analytical solution of the corresponding master equation of such a system. Further, even if the total Hamiltonian 
of an open system-a system plus an environment, is time-independent, the master equation of the reduced density 
matrix of the investigated system still becomes nonautonomous under the non-Markovian dynamics ^(J- Therefore, 
quantum master equation of the reduced density matrix, in general, should be nonautonomous. 

In the previous works [Til IT^ , , we have proposed and developed an algebraic method to treat the sandwich 
terms in the Liouville operator for the nonequilibrium quantum process. This method is just a generalization of the 
algebraic dynamical method |l4j from quantum mechanical systems to quantum statistical ones. According to the 
characteristic of the sandwich terms in the Liouville operator, the left and right representations 0] of the relevant 
algebra have been introduced and the corresponding composite algebra has been constructed. As a result, the master 
equation has been converted into a Schrodinger-like equation and the problems can be solved exactly in the framework 
of algebraic dynamics. This method is very effective to treat the nonautonomous quantum system which contains 
time-dependent parameters for control of the system. 

In this paper, using the algebraic dynamical method, we shall solve the problem of two-level atom interacting with 
the squeezed vacuum field reservoir and investigate what effect on the system can be induced by the squeezed vacuum 
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field reservoir. As is well known, when the environment is ordinary thermal equilibrium fields, the model, either be 
autonomous or be nonautonomous, has been investigated very well H,0|. When the environment is the squeezed 
vacuum field, some work have been done on this field Different to Ref. 0, where two-level atoms with a 

time-dependent external classical driven field in the squeezed vacuum field reservoir were investigated, in this paper, 
we will concentrate on the behavior of a two-level atom in a time-dependent squeezed vacuum field reservoir. With the 
help of the algebraic dynamical method, the su(2)®su(2) dynamical symmetry of the nonautonomous master equation 
for the two-level atom is found. The analytical solutions, both the steady solutions and the time-dependent solutions 
of the system are obtained exactly, and the decay property of the atom is investigated. For the asymptotic behavior 
of the system, it is proven that any time-dependent solution of the system approaches its unique steady equilibrium 
solution. Based on the analytical solution, the different response behaviors of the system to the time-dependent 
squeezing parameter r, linear and nonlinear, are investigated numerically. We noticed that in the nonlinear response 
regime, the expectation of one of the squeezed components of the system, a v , is much enhanced, and the squeezed 
property of the system, i.e., the asymmetry of the decay of the expectation of a x and a y , can also be manifested even 
when the initial values of <7j (i — x, y) are zero. 

The paper is organized as follows. In section II, the model Hamiltonian of the system is presented and the master 
equation for the reduced matrix of the atom is deduced. In section III, the dynamical su(2)k © su(2)j algebraic 
structure of the Liouville operator of the master equation is found and the dynamical symmetry of the system 
is thus exposed. Section IV is devoted to obtain the analytical non-equilibrium solution of the nonautonomous 
master equation, and the approach to the unique steady equilibrium solution asymptotically is proved. In section V, 
numerical results are presented for illustration of non-equilibrium process of some physical quantities. Discussions 
and conclusions are given in the final section. 



II. A TWO-LEVEL ATOM IN THE SQUEEZED VACUUM FIELD RESERVOIR 

The Hamiltonian of a two-level atom interacting with a squeezed vacuum field reservoir is 

H = ^-haj <J z + h^2oj k {ala k + 7:) + h'^g k {(j + a k + h.c), 
1 k 1 k 

Since the squeezed vacuum field reservoir is kept fixed, the total density operator of the system can be written as 

PT (t)=p(t)®l[s k (0\o k )(o k \st(0, 

k 

where S k (£) is the squeezing operator and reads as 

S k (£) = exp(f a ko+k a ko - k - ^4 +fe4 -fe) 

with lo = cfco and £ = rexp(z#), r being the squeezing parameter and 9 being the reference phase for the squeezed 
field. Then the master equation for the reduced density matrix p (t) of the atom interacting with the squeezed vacuum 
reservoir can be obtained with the standard Markovian approximation 0, 

p = -(N+l) (2er_pcr + - (7 + cr_p- pcr+(T-) 
+ — N (2c +( o<7_ — er_<7 + /9 — pcr_(T + ) 

— 7Mcr_ j oa'_ — 7M*(7 +( 0(T + , (1) 

where {a\a k t) = N5 k k> = smh 2 (r)S kk > , {a\a\,) = -M8 k >^k a -k = - cosh (r) sinh (r) exp (-i9) S k >,2k -k, and N is very 
large. Here we have assumed that the vacuum field reservoir is the ideal squeezed one, that is, M 2 = N(N + 1). 
When N — > n —< ^ ><< 1, M — > 0, the master equation JTJ reduces to the familiar form of the master 

ox P(FiT) + 1 

equation 13] describing an two-level atom coupling to the ordinary thermal equilibrium radiation field. The last 
two terms of the above equation exhibit the phase-sensitive nature of the investigated system. The autonomous case 
of this equation can be found in |16(, while in this paper we concentrate on the nonautonomous case where all the 
parameters M, N and, 7 of Eq. are time-dependent, which allow for adjustment of the squeezing parameter r, 
the reference phase 8, and the coupling parameter g k during the time. 
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III. DYNAMICAL SYMMETRY OF THE MASTER EQUATION 

First we will explore the algebraic structure of Eq. (JTJ. Based on the left and right representations of certain 
algebra we can S e ^ the right and left representations of the algebra su (2) = {er+, er_, a z } 

su(2) r : [a r z ,a r ± ]=±2a r ± , [* r + ,a r _] = <r r z , 

su(2) t : [a l z ,a l ± ]=T2a l ± , [a l + ,a l _] = -c/ z . (2) 

It is evident that su(2) r is isomorphic to the su(2), while su(2)i anti-isomorphic to the su(2). This is because su(2) r 
operates towards the right on the bra space | ); but on the other hand, su(2)i operates on the ket space ( |. Just for 
the reason that su(2) r and su(2)i operate on different spaces(the dual bra and ket spaces), they commute with each 
other, i.e. 

[*u(2) r ,m*(2),] = 0. (3) 
From these basic representations of su (2) we can constitute two composite algebras 

su{2) : {J = &±Zl,J+=o%oL,J-=o>Lo 1 + }, 
su(2) : {K = ^^,K + = a r + a l + ,K_=a r _a l _}. 

According to Eqs. J5J and Eq. J3J) it's straightforward to check the following two su(2) commutation relations 

[J 0) J±] = ±2J±, [J+,J_] = J , 
[K 0) K ± ] = ±2K±, [K + ,K^]=K Q , 
[J h K 3 ] = (i,j = 0,±). 

These two su(2) generators have the action on the bases of von-Neumann space, which span the basis of the atomic 
density matrix 
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where s(s ) = ±1. 

By virtue of the above algebras, the nonautonomous master equation JQ) can be rewritten as a linear combination 
of these generators 

p = 7 (t) {[N(t) + 1] J + N(t)J + 1 J - M{t)K. M*(t)K + 2N( ^ + 1 } P = Tp. (5) 

which implies that Eq. possesses an su (2) 7 ® su (2) K dynamical symmetry. Thus it is integrable and can be solved 
analytically according to algebraic dynamics Moreover, the master equation Q is converted into a Schrodinger- 
like equation (0 , where the rate operator T plays the role of the Hamiltonian and the reduced matrix plays the role 
of the wavefunction. 

It is noted that Eq. (JSJ) is a time-dependent generalization of Eq. Q and it is still under the Markovian approxi- 
mation. That implies a basic assumption: the time dependence of the parameters of the master equation don't alter 
the base of the Markovian approximation-the weak coupling assumption of the system and the reservoir. 
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IV. EXACT SOLUTION TO THE MASTER EQUATION IN THE NONAUTONOMOUS CASE 

A. Steady solution of the master equation 

To better understand the time-dependent solution of the master equation and its decay behavior, we first consider 
the long-time case (7 (t) — ► 7, N(t) — > N,M(t) — > M) and the eigen equation problem of the master equation JSJ. 
From Eq. (J5J we can get 

Tp = (3p (6) 

Introducing the similarity transformation 

U — e ot+J+ e a — T- e ri+K +e ri-K- 

and under the transformation parameter conditions 

(N + l)a\ + a+ - N = 0, 

(A r + l)(l + 2a + a_)+a_ = 0, (7) 

Mr)l - M* = 0, 

l + 2r ?+ 7 ? _ = 0, (8) 

we can transform Eq. © to the diagonal form as 

Tp = f3p, 

p = u-'p, 

f = U- 1 TU = -j{[(N + l)a + + ~]J -Mr)+K + ??^}. (9) 
The eigensolutions of Eq. © are 

p( s ,s') = e a + J +e a - J -e ri + K +e^- K -\s){s'\. (10) 

It is interesting to note that Eqs. Q and Eqs. (JSJ both have two sets of solutions. The two combinations of 
the solutions (a+ = -1, a_ = (TV + 1)/(2N + 1), 77+ = ±e ie , and ?7_ = Te~ ie ) and (a+ = N/(2N + 1), a_ = 
— (AT + l)/(2N + 1), ??+ = ±e zB , and ?y_ = =Fe - * 6 ') both contain the same zero-mode solution (the unique steady 
equilibrium solution) and they are equivalent physical solutions since all their eigen values 0{s, s') are non-positive 
(this can be proved by the series expansion M = y / N(N + 1) w N(l + ^ — §^7) and by the fact N ± M + | > 0). 

Noticing that J± yield non-zero results only if they act on diagonal elements and K± do so if they act on non- 
diagonal elements, we see that the zero mode steady solution of the system has the same form (but with different N) 
as that in Ref. [13J where the reservoir is a thermal equilibrium field, 

N + l , N , w 

^ = ^ TT l-i)(-i| + ^TTl + i)(+i| ! (ii) 

It is noted that the rate operator T is non-Hermitian, i.e. 1^ 7^ T, which is evident from jt = J-, J_ = J+, 
Jq = Jo, K' + = K_, and K_ = K + . Just because of this non-Hermiticity, the eigenvectors of V and constitute a 
bi-orthogonal basis [l8|. By introducing a similarity transformation U' — (U^ 1 )^ and umder the same conditions as 
Eqs. J71IHJ), the operator can be diagonalized. Then the eigensolutions of are given by 

rtp(a, a ') = p(s,s')p(s,s') 

p{s,s') = e- , '+ K -e-^ K +e- a + J ~e- a - J +\s)(s'\. (12) 



One can easily to check that the eigenvectors p(s, s') and p{s, s') form a biorthogonal set. A similar discussion can be 
found in Ref. 0. 
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To compare with the work of Ref. pH. I20I l2l| . we consider the left eigensolutions of T. The transformed form of 
the eigen equation is 

Eq. © shows that the operator f = U~ X TU is self-adjoint and, thus, p'T = WT^l/- 1 )^' = [(J/- 1 )*]" 1 ^^- 1 )^'. 
This implies that T' has the same eigensolutions as T when they act from the left and right on p, respectively. In the 
phrase of Ref. Jig , p(s, s') is the right-hand eigenvectors of the operator V and p(s, s') is the right-hand eigenvectors of 
the operator T'. Furthermore the right-hand eigenvectors of the operator coincide with the left-hand eigenvectors 
of r. And the two sets of eigenvectors constitute a biorthogonal set. From this biorthogonal set one can construct a 
generalized entropy functional varying monotonically with time, which named as the Lyapunov functional and have 
a crucial rules in control theory 20] . 



B. Time-dependent solution of the master equation 

Now we turn back to the nonautonomous case and study the time-dependent solution of Eq. JSJ. With the gauge 
transformation 

U g (t) = e Q +(*) J + e Q -W J - e v+(t)K+ e v-(t)K- ^ 

and the gauge parameters satisfying 



<ia+ (t) 

dt 
da- (t) 

di 
dy+ (t) 

dt 

dy- (t) 

dt 



-l(t)[N(t) + l]a\(t) - 7 (t)a+(t) + 
y(t)[N(t) + 1][1 + 2a+(t)a-(t)} + 7 (t)a_(t), 
7 (t)[M(t)rjl (t)-M*(t)}, 

- 1 (t)M(t)[l + 2 V+ (t)^(t)}. (13) 



where the initial conditions of the gauge transformation is U g (0) = 1, the operator T can be transformed to the 
diagonal form 

f(t) = U^{t)TV g {t)-Uj l (t)U B tt) 

= -j(t){[(N (t) + l)a+ (t) + \]Jo~ M (t) v+ (t) K + 2N ^ + 1 } 

If the initial state of the system is p (0) = ^ s , \ s . s > \s)(s'\, then the time-dependent solution of the nonautonomous 
master Eq. JSJ is 

P (t) = ^A s , s ^ s (i)e^-^ T ){[^W+ 1 ) a +W+i]^- M ^+M^+ M ^> dT | S )( S / | 

8,8* 

= {Ai,i/i,i(t)[l + a+ (t) a_(t)] + A_i,_i/_i J -i(t)a+ (t)}| + 1)(+1| 
+[Ai,i/i,i(*)a_(*) + A_i,_i/_i,_i(t)]| - 
+{Ai,_ 1 /i,_i(t)[l + y + (t) r,-(t)] + A_ 1 , 1 /_ 1)1 (t)7 7+ (t)}| + 

+ [A ll _i/i,_i(t)»/-(*)+A_ 1) i/_i, 1 (t)]|-l)(+l|, (14) 

where we have defined f„, s >(t) = eti -7(r){[(JV(r)+i) a+ (r)+i]^~M(r)^ + (r)^+ 2JV( ; )+1 }dr^ 

One can check that the solution of Eq. i|14|) recovers the results of the autonomous system with time-independent 
parameters. For the autonomous case, the explicit solutions of the gauge parameters a± and i]± can be obtained from 
Eqs. O 

1 _ - 7 (2JV+l)i 

a + (t) 



N+l , p - 7 (2JV+l)t ' 
N ' 



«-(*) 



(N + 1)N[^ + e -7(2Af+i)t][i _ e -7(2iV+l)i 

(27V + l)2 e -7(2W+l)t 
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1 - e 2 ~< Mt 



1 + e 2 -r Mt ' 
(1 - e 2 T Mt )(l + e 2 i m ) 



4 e 2 7 Mt ' 

where M — M* is assumed. Then under the initial condition of the system 

p(0) = | M | 2 |1)(1| + W\ 2 \ 1}{-1| + /*f*|l><-l| + M*H - 1)<1|, (15) 
the expectation values of a x , a y , and a z are 

(o») = 7(M^-^)^ (JV - M+I)t , 



2[| / ,| 2 (jV + l)-|H 2 iV]e-^ jv+1 ) t -l 
(*.) - ■ (16) 

So the decoherent characteristic times of (cr x ), (<Ty), an d (<j z ) are 7 ^ N+ \ I+ 1^ , 7 ( Ar _^, /+ i i an d 7 ( 2 at+i) ' res P ec tively, 
which is same as the results of Ref. 0|. The decoherent characteristic time of (a y ) is much larger than that of 
(a x ). From the above results we see that the radiative properties of the atom depend sensitively on the state of the 
environment to which it is coupled. The squeezed vacuum reservoir imposes its phase information on the atom and 
makes a highly asymmetric decay property of the atom. This is quite different from thermal vacuum reservoir which 
results in the same decay rate (where M = 0, N — ► n) for (a y ) and (cr x ). 

So we can see that the time-dependent solution of Eq. I|14|) reproduces all the results of the autonomous case. In 
the following we shall concentrate on the nonautonomous case and investigate the time evolution behaviors of the 
system. 



C. The asymptotical behavior of the solution 



In this section we shall examine the asymptotic behavior of the time-dependent analytical solution (|14|) of the 
nonautonomous master equation and prove that it approaches to the equilibrium steady solution To this end 

we shall first investigate the asymptotic behavior of the gauge parameters which are solutions of Eqs. Ijl3|) from 
which one notices that drj + (t)/dt < 0(> 0) if -1 < r) + (t) < (if r]+(t) < -1 and r) + (t) > 0). With the initial 
condition r/ + (0) = 0, we see that r] + (t) approaches the value 77+ (00) — — 1 asymptotically. To get the asymptotic 
solution of ?7_ (t), we define y(t) = ry_ (t) eJo 2 7 {r)M(r)r, + (r)dr _ rj _ ^ e ,Q P {r)dr_ The time differential of y ^ is given 

by &exp J* p(t)cIt where b = f]-{t) + r]-(t)p(t). Since b — ► 7 which is bounded and p(t) is negative for large t, 
the differential y (t) tends to zero, hence y (t) tends towards a constant. With the same procedure we can get the 
asymptotic solutions of a± (t) . Summarizing the above discussion, we have 

N 



N+l 

a_(t) x e-tf ■y(r)[N(T)+i][a + ( T )+i]dr^ oo = cmst 

V+ (*) It^oo = -1 
^(i^o^MMM^Mtfri^ = C(mst 

Although, for the steady case, two sets of the parameter solutions a± and rj± lead to two equivalent transformations 
which diagonalize the rate operator and generate the same physical solution, as discussed above; for the time-dependent 
case, the properties of the two sets of steady solutions of a± and r]± are not on equal footing. It is found that only one 
set of the solutions can be reached by the time-dependent solution asymptotically. Using these asymptotic relations 
we have the asymptotic results of the time-dependent solution as follows 

P++(*)l*-=o = e-/o 7(r)[JV(r) + l][a + (r) + l]«ir [ | + + Q _ (t)( | _ + 1> <+l|)] 

— > const, x p 

P—(t)\ t ^oo - e -/o*7(r)[-(Bo(r)+l)a + (r)+fio(r)]«Ir[| _ 1) (_ 1 1 + a+ ( t ) | + + _^ const . X p 
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/* 7 (r)[M(r)^ + (r)+Ar(r)+|]dr e / '2 7 (r)M(r) J)+ (r)dr { | + !)(_!| +r; _ (;) ( t ) | + 1> <_ 1 1 + | _ 1>( + 1|]} 
/ ' 7 (r)[M(r) + ;V(r)+i]dT| + + g - J* 7 ( T ) [JV (r)-M(r)+i]dr x confit[ | _ + _ | + 

J ' 7 (r)[M(r)^ + (r)+Ar(r)+i]dr^ + ^ | + + | _ 



The above results indicate that any time-dependent solution of the master equation in the nonautonomous case 
asymptotically approaches the unique steady equilibrium solution, irrespective of their initial conditions. 

V. NUMERICAL RESULTS 

To study the squeezed properties of the system in the nonautonomous case, we shall start from the time-dependent 
analytical solution Eq. (|14|l and calculate the expectation values of the two-level atomic operators. The expectation 
values oi <Ti (i = x,y, z) are 

, <, + V+(t)y-(t) + 77- (t)] + + y+(t)] 
K a x) = 2 ' 

, \ nv*f 1 ,-i[l + r,+ (t)T,-(t) -y-(t)} +//V/_ 1 ,i[7 ?+ (t) - 1] 
K> = 2l ' 

(<Tz) = |H 2 /n[l + a+(i)«-W-"-W] + H 2 /-L-iKW-l], (17) 

where the initial state is same as Eq. (|15|) . Given any time-dependent parameters "f(t), r(t), and 9{t), we can get 
the gauge transformation parameters a±{t) and rj± from Eqs. I|13|) numerically. Then substitute a±(t) and r]± to 
Eqs. H17[l. we can obtain the time-dependent behaviors of the expectation values varying according to the system 
parameters. In the following we will study the variation of the time-dependent behavior of (<7j) with the squeezing 
parameter r(t). 

Without loss generalization, to explore the time-dependent response of the system to the squeezing parameters we 
choose time-dependent parameters as r = cie _C2 *, 7 = 1, and 9 = 0. Fig. 1 and Fig. 2 show the time-dependent 
behaviors of {at) with ci = 0.1 and C2 = 0.1 under the initial condition /1 = \/0.2e 1 / 3 , v = \/0.8e 27rl and /1 = V0.2, 
v = V0.8 respectively. Fig. 3 and Fig. 4 correspond to C2 = 0.3. The behavior of (cr^) varying to r(t) with parameter 
C2 = 0.6 were plotted on Fig. 5 and Fig. 6. From these figures we see that (<r z ) decays to its steady value — 1; 
(a x ) and (a y ) decay to their steady value 0, but the decay time of (a y ) is much larger than that of {(J x ), which 
regenerate the results of the autonomous case. From Fig. 2 we notice that (<r y ) = 0, so the squeezed property of 
the squeezed vacuum field to the two-level atom is covered by the initial value of {o y ). For the small squeezing 
parameter, the time-dependent response of the system with respect to the squeezing parameter is nearly linear and 
much like the autonomous case. However, with the increase of squeezing parameter, the response is no longer linear 
and the decay time of (<7a;)'s becomes shorter and shorter, while that of (cr^'s becomes longer and longer. A more 
interesting phenomenon appearing during the increase of the squeezing parameter is that the value of (a y ) can be 
much enhanced. And the resulting effect of the squeezed vacuum field reservoir on the two- level atom, i.e., the 
large asymmetric behaviors of decay between (a x ) and {o y ), can also be presented even when some of them are zero 
initially, which is not presented in the autonomous case (this can also be seen from the second expression of Eqs. 
mot h Moreover, the time-dependent behavior of (a y ) is completely different from that of r(t), and the response of 
system to the squeezing parameter thus exhibits nonlinearity. 

From the analysis above we see that the squeezed vacuum field reservoir acts as a squeezing mold and imprint 
its squeezing information to the two-level atom, which makes the asymmetry of the two decay behaviors. With a 
time-dependent squeezing parameter this action can be more enhanced, which widen the results of Ref. [HI and 
provide a potential way to get more stable atomic variable. The similar time-dependent system can also be found in 
Ref. 22], where a technique of stimulated Raman adiabatic passage was used to get a complete population transfer 
between two atomic states. There the two adiabatic varying lasers act as a role of coherent trapping in the three-level 
atomic system and the technique relies on the coherent trapping state. 

VI. SUMMARY AND OUTLOOK 

We have investigated the nonautonomous master equation of the two-level atom interacting with the squeezed 
vacuum reservoir for the first time. The new results of our paper compared to the previous ones are: (1) the 



p+-(t)\ t ^oo = e 

— > 

P— l-(*)U— 00 = e 

— > 
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dynamical symmetry of the nonautonomous master equation has been established and the master equation has been 
converted into a Schrddinger-like equation. Based on the dynamical symmetry we have obtained the most general 
non-equilibrium solutions of the nonautonomous master equation analytically; and as a special case, our solutions 
recover all the results of the autonomous case. (2) We have proven that any time-dependent and non-equilibrium 
solution of the nonautonomous master equation approaches its unique steady equilibrium solution asymptotically. 
(3) The non-equilibrium solution of the nonautonomous master equation is analysed numerically and the different 
response behaviors of the system to the squeezing parameter is exposed, and nonlinear effect is exhibited. 

The works of dissipative system combined with the research of control theory were extended in some publications 
|20l l2l| , how to connect our approach with this investigate is a interesting and challenge problem. 
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1: The time-dependence of with parameters r = O.le 01t , ji = V0.2e 7rl,/3 , and ^ 
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FIG. 2: The time-dependence of (cri) with parameters r — O.le 01t , ^ = \/0.2, and i/ 



A 

V 



0.4 -, 

0.2 

0.0 
-0.2 -| 
-0.4 
-0.6- 
-0.8- 




-1.0 



(«**>) 
(«*>) 
(<">) 



10 



20 



FIG. 3: The time-dependence of with parameters r = 0.3e 01t , ji = V0.2e 7rl,/3 , and 



A 

b 

V 



0.4 -, 



0.2 



0.0 



-0.2- 



-0.4 



-0.6- 



-0.8- 



-1.0 



(<*>) 
(<^ y >) 
(«*>) 



10 



20 



FIG. 4: The time-dependence of (cr») with parameters r = 0.3e ^ = \/0.2, and 



